On the Integrability of the 72-Centre Problem 

Andreas Knauf* Iskander A. Taimanov^ 



Abstract 

It is known that for n > 3 centres and positive energies the n- 
centre problem of celestial mechanics leads to a flow with a strange 
repellor and positive topological entropy. 

Here we consider the energies above some threshold and show: 
Whereas for arbitrary g > 1 independent integrals of Gevrey class g 
exist, no real-analytic (that is, Gevrey class 1) independent integral 
exists. 



1 Introduction 



In |BTj the existence of a smoothly integrable geodesic flow on a compact 
manifold with positive topological entropy was established. Positivity of 
topological entropy is seen as an indication of complex dynamics, whereas 
integrability of a Hamiltonian flow is a metaphor for its simplicity. So coex- 
istence of these two properties may not have been expected. In fact, as we 
show here, such a coexistence takes place in natural physical problems. 

To be more specific, in this note we consider the n-centre problem of 
celestial mechanics in d = 2 and 3 dimensions. We denote by Sk G Mf^, Zk G 
M \ {0} the location resp. strength of the fc-th centre, assuming Sk 7^ si for 
k ^ I. Then the Hamiltonian function 

H ■.T*M -^R , H{p,q) = \f + V{q), 
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with potential 



V:M^R , V{q) = -J2j^^ (1) 

k=i ^'^ll 

on the cotangent bundle T*M of configuration space 

M :=M'^\{si,...,4} 

generates a - in general incomplete - flow. 

Denoting by a) := J2i=i ^ ^Pi \t*m restricted canonical symplectic 
form, there exists a unique smooth extension 

{P,u!,H) of the hamiltonian system {T*M,(2},H) (2) 



such that the flow $:MxP^P of if is complete (see jKnj . Thm. 5.1). 
Concerning integrability of the flow, the following is known: 

• For n = 1 this Hamiltonian system is integrable, with angular momen- 
tum 



L: P ^R-^ , L(x) 



{q- si) xp , x = {p,q), q^ Si 
, otherwise 



for dimension d = 3 being a real analytic constant of motion. For 
Zi > this is called the Kepler problem. 

For n = 2 centres one introduces elliptic prolate coordinates to analyt- 
ically integrate the flow $, see e.g. |Arj or |Thj . 



For n > 3 centres and d = 2 Bolotin showed in jBoj the nonexistence of 
an analytic integral of the motion which is non-constant on an energy 
shell H~^{E), E > 0, see also the discussion in Fomenko |Foj . 

For d = 3 and a collinear conflguration of centres the angular momen- 
tum w.r.t. that axis is an additional constant of the motion, indepen- 
dent of the number n of centres. However, for = 3 it was proved 
by the flrst author |Knj for sufficiently large energies E > Et^ and by 
Bolotin and Negrini |BN1 IBN2j for nonnegative energies E > that 
the topological entropy of the fiow, restricted to the set 6^; of bounded 
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orbits on H~^{E) is positive (and htop{E) = for Be = 0)- Further- 
more htop{E) is zero for n = 1 and 2, and htop{E) > if n > 3 and all 
centres being attracting or not more than two Sk being on a line (for 
collinear configurations with Zi, . . . , Zn < one has htop{E) = for 
E > 0). 

In the present paper it is proved that 

• for d = 2, attracting centres (Zk > Oj and E > 0, 

• for d = 3, arbitrary Zk ^ 0, non-coUinear configurations of centres and 
E > Etii where the threshold energy level Eth is determined by the 
data Zk, Sk, k = 1, . . . ,n, 

the n-centre problem restricted onto the energy level E admits d—1 indepen- 
dent integrals of motion which are smooth (and moreover are of the Gevrey 
class g for any g > 1) (see Thm. [^. Therefore the restricted problem is 
smoothly integrable. 

On the other hand if the affine span of the (non- collinear) centres equals 
M.^ then the restricted problem is not integrable in the real-analytic sense, 
that is Gevrey class 1 (see Thm. [2J). 

The article is based on the analysis of the n-centre problem given in the 
paper by the first author |Knj where, in particular, it is shown that for high 
energies there are maps which relate scattering orbits to their asymptotics 
given by scattering orbits of the Kepler problem which is integrable (the 
necessary facts extracted from |Knj are exposed in sections 2-3 of the present 
paper). Such a relation of the n-centre problem with an integrable problem 
leads to the smooth integrability of the n-centre problem which is proved here. 
Here we use a trick for constructing smooth first integrals from discontinuous 
preserved quantities similar to others used in |But IBTj 

We would like to notice that in |Boj the nonexistence of an additional 
analytic integral of motion for the two-dimensional problem was derived by 
Bolotin from Kozlov's theorem |Koj which reads that the geodesic flow of 
(real)-analytic metric on a closed oriented surface of genus g > I does not 
admit an additional analytic first integral. 

Although one can expect that the integrability of a problem by analytic 
integrals of motion implies vanishing of the topological entropy it is not 
proved until now. Therefore the results from IKnj do not imply the 
nonexistence of a complete family of analytic first integrals. We prove that by 
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an analysis of the set formed by bounded orbits which supports the restricted 
flow with positive topological entropy. 



2 Known Smoothness Results 

There are three basic types of motion: bounded, scattering and trapped, 
corresponding to the disjoint subsets b,s,t C P with 

:= {x e P \ g (±M+, x) is bounded } , 6 := 6+ n 6", 

s± := {x e P I X ^ 6=^ and H{x) > 0} , s := s+ fl s" 

and t := s^As~. The orbits in s go to spatial infinity in both time directions, 
but because of the long range character of the effective potential of strength 
=1 Zk, one describes these limits by comparison with (regularized) 
Kepler flow $^ : Poo — ^ Poo generated by the extension of 

H^:T*iR'\{0})^R , H^(p, q) -.= Ip^ - (3) 

Identifying the two phase spaces P and Pqo outside a region projecting to a 
ball in configuration space which contains all singularities, the M0ller trans- 
formations 

n± : Poo,+ := {xeP^l H^{x) > 0} , := lim $-* o Id o (4) 

exist as pointwise limits, and are measure-preserving diffeomorphisms (Thm. 
6,3 and 6.5 of |Knj ) . Similarly the asymptotic limits of the momentum 

-.s^ ^R'^ , p^ixo) := lim po $*(xo) (5) 

are smooth. Finally we define time delay r : s — > M of a scattering state 
X E s by 

t{x) := lim ! (a{R) o - (6) 

lia^R) o o ntix) + a^{R) o $^ o fi;(a;)))cit, 

where cr(P) : P ^ {0, 1} and croo(P) : Poo {0, 1} are the characteristic 
functions a{R){x) := 9{R — ||^(a^)||) and similarly for croo(P). That asymp- 
totic difference between the time spent by the orbit and its Kepler limits 
inside a ball of large radius diverges near but. However r is smooth, as the 
M0ller transformations are. 
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3 Analyticity Properties 

Proposition 1 For a potential V of the form (QJj the following maps are 
real- analytic: 

1. the flow $ : R X P -> P 

2. the M0ller transformations fi^ : Poo,+ and asymptotic momenta 



3. the time delay r : s — M. 
Proof. 

1) We first indicate the definition of phase space P in order to show that for 
the potential (0) the smooth extension (0) actually works in the real-analytic 
category (in |Knj more general non-analytic potentials V were considered). 
We assume li = 3, the case oi d = 2 dimensions following by restriction. 
For small e > in the phase space neighbourhood 



of the /th centre the following real-analytic coordinates are used: 
• The restriction of the Hamiltonian function, which splits into 






H(p,q)=Hi{p,q) + Wi{q) 



with Hi{p,q) := ^p^ 



q- si 




and 




• The angular momentum 



Li : f/f ^ M' 
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Li{p,q) := (g - si) x p. 



(9) 



relative to the position s/. 
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The time Ti : f/f ^ R after which the Kepler orbit generated by Hi is 
in its pericentre w.r.t. Si 

Ti{p,q):= I ■ sign((g- g;) -p). 

Jr^Um ^2r'^Hi{p, q) + 2Z;r - L'j{p, q) 

(10) 

The neighbourhood Uf C T*M is defined in a way which makes the 
pericentre unique. Its distance from si equals 

TnAMq) ■■= { . 2H,(P,9) '"^'^^ . (11) 

Lf{p,q)/2Zi ,Hi = 



The Runge-Lenz vector relative to si is given by 



Fr.U^^R' , Fi{p,q):=pxLi{p,q)-Z, 



Since \Fi\'^ > Zf /A > 0, the pericentral direction 

Fr.U^^S' , Fr.= Fi/\Fi\ (12) 

is well-defined. 

As ■ = 0, we get a real- analytic diffeomorphism onto the image 

y:U[-^ T*(R X S') \ 0, {p, q) ^ {fi, U; Hi, Fi) 

onto a punctured neighbourhood of the zero section of the cotangent bundle. 
The fiow generated by Hi is linearized in these coordinates, and extended on 
Uf ■= JJf U (M X ^2) by 

3^:=(r„L,;ff„F,):C/f^r-(MxS=),3^W:.nO.O;ft,f) (/.f).MxS= 



y{x) otherwise 

to a full neighbourhood of the zero section. 

In a final step the hamiltonian fiow $ : M x P ^ P generated the con- 
tinuous extension if : P ^ M of the Hamiltonian function H : T*M M 
is linearized by slightly changing the coordinates y using the (real-analytic 
but incomplete) fiow generated by H. See |Knj for details. 
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In summary, we obtain a real-analytic extension Q of the Hamiltonian 
system {T*M ,(l!,H). Thus the Hamiltonian vector field Xh (defined by 
ixff^ = dH) is real-analytic, too, and thus the fiow $ : M x P ^ P is known 
to be real- analytic (see, e.g. |Hoj ) . proving assertion 1. 

2) This, however is insufficient to show real-analyticity of the M0ller trans- 
formations and asymptotic momenta. It is known that even for smooth 
potentials these maps may be very non-smooth, see jSij. 

There exists an energy- dependent virial radius Pvir > 0, with 

^ {q.V) > f > if M\ > R^iriE) and E := H{p,q). (13) 
If we assume go '■= WoW > -Rvir(-E) and (goiPo) > 0, then 

\\q{t)\\>(lo- {t)x for alH > 0, with X:=^/E/2/qo (14) 

and 

{t), := Vl + (At)2 , (t):=(t)i. 

In particular a trajectory leaving the ball of radius Pvir(-£') cannot reenter 
this ball in the future but must go to spatial infinity. 

Analyticity estimates for the trajectory are then derived from the integral 
equation with initial conditions Xq = {po,qo) 

q{t, Xq) =qQ + tpo- / VF(g(r, a;o)) dr ds, 
Jo Jo 

using the decay property of the potential 

/ ^ \ 1/31+1 



outside the virial radius, and (fT^ . 
Setting 

\\w\\\ := sup 



t>o {t)x 

we obtain for 7 := {a,/3), := (9" 9^ 



||92„g(-,xo)||A < a!go''^"'''^'^^i5-^'"'-'+''^i'^ (15) 
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This allows us to perform the (locally uniform w.r.t. Xq) time limit in 

dl{p{t,xo)-po) = - (16) 

D^VV{q{T,Xo))[dt;W{ T,xo), . . . ,d2g q{T,xo)]dT. 



N=l ^(i) + ,..+^{iV)^^ 
|7(')|>0 



with g := I7I > 1, and to conclude that p"*" is real-analytic at xq. We can 
substitute the assumption xq G for the stronger assumptions go := Ikoll > 
R^i^{E), {qo,po) > 0, as initial conditions meeting the first lead to data 
meeting the second one after some t. The same holds for p~ by reversibility 

{p~{po,qo) = -p+(-Po,go))- 

The proof of real-analyticity of the M0ller transforms is based on the 
integral equation 



POO 

It J 









Qirixo] 



— V{q{T, Xq)) drds 



for r(t; a;o) := q{t,Xo) — Q{t;xo) with Kepler hyperbola yP{t; Xo),Q{t; Xq) 

^^(Xo). Inspecting the proof of smoothness for the M0ller transforms in 
[Kn , one sees that the estimates for dj^^r can be dominated by |7|!Cl'''l. 

3) For a scattering state x E s time delay equals 

t{x) = ilim /"cr+o$*(x)- (aoo(i^)o^]+-aoo(i?)o^];)($*(x))rf^ 

+ ilim /"(T-o$*(x)- ((Too(i?)ofi;-aoo(i?)of]+)(<l>*(x))c/t, 

with cr^(p, q) := 9{±q-p), see |KKj and Qf := (i7^)~^. Using the intertwining 
property Qf o = <|)^ o this can be written entirely in terms of M0ller 
transformations and the Kepler flow: 



rix 



i hm / a+ o ^l{x) ■ {a^iR) o o ntix) - a^iR) o o n;{x)) dt 

- 1 hm / a- o $*^(x) ■ (a^R) o o fi+(x) - a^R) o $^ o fi;(x)) rft. 
With Assertion 2 and the analog of formula ()10|) this implies analyticity of 

T. □ 
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4 Gevrey Integrals of Motion 



We now show the existence of independent constants of motion for all energies 
E > Eth- In |Knj many estimates are shown to hold true above a threshold 
energy -Eth > 0: 

In |Knj many estimates are shown to hold true above a threshold energy 
Eth > 0: 

• For d = 2 and attracting centres {Zj. > 0) i?th = 0. 

• For d = 3, arbitrary and non-coUinear configurations of the 
Sk G the existence of such a threshold is proven. 

In particular the set Be of bounded orbits of energy E > E^y^ is shown to be 
of measure zero (and has a Cantor set structure for n > 3). 

According to the standard definition (see, e.g. jAMj . Def. 5.2.20) functions 
/o, . . . , /fc : P — M on a symplectic manifold (P, uj) are called independent if 
the set of singular points of F := /o x . . . x : P ^ M'^ has measure zero. 

To simplify discussion, we set 

P:=H-\[E,,E2]) 

for an arbitrary energy interval, Eth < Ei < E2 < oo, and /q = /o := H\p. 
Then for a parameter > 1 we define for A; = 1, 2 



(17) 



where C{g) = with C2 to be defined below. Putting things together we 
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get a map : P R^. 

The notation is chosen so that F^ belongs to the Gevrey class of index g. 

Now we collect some information about Gevrey functions, which were 
introduced in 



Definition 1 For g > I and an open set C a function / G C°°(fi) is 
called of Gevrey class g if for every compact K C Q there exist Ak, Ck > 
with 

max,g;,|9°/(x)| < AkC^k^ {a\y (a E N^). 

Here |a| = ai + . . . + q;„ and a! = ai! ■ . . . ■ aj. if a = (ai, . . . , a„). 
The vector space of these functions is denoted by Gg{Q). 
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Then Gi{fl) is the space of real- analytic functions, and Gg>{Q) D Gg{Q) 
for g' > g. Gg{Q) is stable w.r.t. partial derivatives, compositions, and the 
implicit function theorem holds within the class. For g > 1 Gevrey partitions 
of unity exist. 

Example 1 For g > 1 

._ / exp(-x-i/{^'-i)) , X > 
^ 1 , X < 0, 

is a function in G'g(M) but not in ^^/(M) for any g' < g, see e.g. jJuj. 

For a real-analytic manifold we can define the spaces of Gevrey functions, as 
the defining family of bounds is preserved by coordinate changes. 

Theorem 1 For all g > 1 the functions Jq, ff, /| are independent, in invo- 
lution and of Gevrey class g. 

Proof. 

In Thm. 12.8 of |Knj it was shown that the set Be of bounded orbits of 
energy E is of Liouville measure zero for all E > Eth- As the set t^; := {x G 
H~^[E) I X G T} of trapped orbits is always of Liouville measure zero, on 
P the set P \ s is of measure zero. On P fl s the functions /f and /| are 
real- analytic, using Proposition Q whereas they are zero on P\s. 

We study their decay near P \ s in order to prove that they are in G^{P). 

All orbits in P \ s enter the interaction zone, so we need only orbits in 
P n s entering the interaction zone. W.l.o.g. we assume Pvir to be constant 
on [El , E2] . Then the region 

i:={xeP \ ||g(x)|| <P,iJ 

projecting to the interaction zone is compact. The restriction of the real- 
analytic flow $:MxP^]Rtoa domain of the form [—e, e] x I thus has 
partial derivatives < CiC^^^al. We conclude, using Corollary [T] 

that for arbitrary t G M such that x) G /, too, 

\d"^t,x)\ < Ciexp(C2|a| (t)) a\. (18) 

Next we analyse time delay (jH)) for orbits entering /. By compactness of / 
there exists a P > Pvir such that all Kepler orbits o Qf{x) | t G M} 
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/ 



Figure 1: Scattering trajectory in configuration space (solid line), with asso- 
ciated Kepler hyperbolae (broken lines) 

through points a; G / enter a configuration space region of radius R. In © 
we only consider radii R > R and, denoting by Tfj{x) the integral ©, we 
split this quantity in 

M^) = 'Tnix) + r^(x) + r°(x). 

Here 

r±(x) := / (a+o$,(x)-(T(i?)(l-a(i?.i,))o$*(x) 

-'^a^iR)o^lontix))dt (19) 

and 

r°(x) := [ (T(i?vir) o <|>*(x) dt. (20) 
Jr 

Whereas the limits r^(x) := lim/j^oo ^^(2;) of (fTIH) meet Gevrey class 1 
estimates uniform on Jfls, t^{x) defined in ()20|) . the time spent by the orbit 
in the interaction zone, is not uniformly bounded on / fl s for n >2 centres, 
see Figure 1. 

We need only consider scattering states in a slightly smaller domain := 
{x G P I ||(?(a;)|| < -Rvir — £}■ By (fT!?jl the configuration space trajectories 
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with initial conditions in n s intersect the boundary of the interaction zone 
with angles uniformly bounded away from zero. Therefore (|18p shows that 
for X G /e n s 

|<9"r°(x)| < Ciexp(C2|a| (r(x))) a!, 

and T = + + satisfies the same kind of estimate (with enlarged 
constants). 

A similar kind of reasoning applies to the asymptotic momenta p^: 
\dyt{x)\ < Ciexp{C2\a\ {t{x))) a\ 

for X G n s. 

IFtoiu. this one can conclude that /f and /f are of Gevrey class g for C (g) 
large enough, see Proposition |21 

The functions and r are $t-invariant. Therefore ff and /| are 
invariant, too. As /q = H \ P generates the flow $t on P, the Poisson 
brackets {f^, ff} and {f^, /|} vanish. 

By definition p^{x) = lim^^-i-ooP ° ^t{x), and the decay estimates (fTH|) 
imply that 

\ lim pi o lim p2 o <l>t(x) i = hm{j9i,j92} o ^t{x) = 0. 

Thus the Poisson bracket {ff, /|} vanishes, too. 

To see that the functions are independent, it suffices to show that the 
(analytic) restriction 

^ \pns 

has a measure zero set of singular points. As the energy range [Ei, E2] under 
consideration consists of regular values of H, we may consider instead the 
independence of ff and /|, restricted to H~^{E) fl s, for E G [Ei, E2]. 

Independence of ff and /|, restricted to an energy surface H'^{E) with 
E G [-El, -E2] follows from their real-analyticity on H~^{E) fl s and the inde- 
pendence of pf and P2, remarking that for {q,p) 00, p^{{p, q)) ^ P- □ 

Proposition 2 Consider real- analytic maps 

f-.^f^w" , g-.n.^nj 

with open sets C M*^, meeting 

\^''f^\<Mfdfa\ , <M,Cf Id (z = l,...,ci,aGN^). 
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Then 

\d"{f og)^\< Mf C|^' a\ with Cfg := Cg{l + CfdMg)d. 
Proof. The chain rule has the form 

n 

k=l Ai,...,Ak 

with n = \a\, {Ai, . . . , A^) running through the fc-partitions of {!,..., n} 
and G Nq is the muhiindex of size la^J = \Ai\ corresponding to Ai C 
{1, . . . , n} in a, say lexicographic, ordering of the n partial derivatives. There- 
fore with A := CfdMg 

n k 

fe=l M,-,Ak '=1 

n k 

= MfC^;Y.^CfdMgfk\ J2 n^-^^' 

fc=l Ai,...,Ak 1=1 

n k 

< MfC^Y^x'ki E ni^'i' 

k=l Ai,...,Ak 1=1 

k=l 

since A;! nf=i l^d' — ^' (fc-i) • As in d dimensions 

|a|!<cil"la! (« e N^), 

the estimate follows. □ 
This estimate is iterated in order to estimate the flow for long times: 

Corollciry 1 Assume that for k = 1, . . . ,t the maps 

admit the uniform estimates |(?"//*^^| < MC'"' a\ with M > 1. 
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Then their iterates T^^^ := /^^^ o T^'^ with T^^^ := Idj^^ are estimated by 

|9"Tf ^1 < Mexp(C|a| {a e N^k = 1, . . . ,d) (21) 

with C := \n{d) + MdC. 

Proof. (j2H) holds for /c = 1, and is assumed to hold for some k < t. Then 
by Proposition 121 



|a"Tf+'^| < M{exp{Ck){l + MdC)dp 



< Mexp((:7(A; + l))l"la!, 

proving (j2H). □ 

Proposition 3 For the choice C{g) := in [Tl\) the functions ff. are of 
Gevrey class g. 

Proof. As \d" {t{x)) \ < Ciexp(C2|a| {r{x))) a!, 

|a| k 

\crexp{C{g){r{x)))\<exp{Cig){r{x)))J2C{g)' Yl H 

k=l Ai,...,Afe 1=1 

< exv{(C(Q) + C,.n)(T(x))] I Cfo)'^'^ fr!)) a\ 



,k=l 



= exp {{Cig) + C^n) {r{x)) ) (1 + C{g)C^Y a\ 
with n := 

The (r) dependent part of |i9"/|| is of the form 



exp (-e^(^)<^>) e^^"<^> < 
as it has its maximum for (r) with 

,C(3)(r) _ Can 



Thus the choice C{g) := leads to the proof. □ 
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5 Nonexistence of Analytic Integrals of Mo- 
tion 



We start with the following simple observation for a hamiltonian system 

Proposition 4 Let 7 C := H~^{E) be a periodic orbit of the Hamilto- 
nian flow generated by the Hamiltonian function H which is isolated on 
the energy surface Sg. Then there is no additional integral of motion which 
is functionally independent of H on 7. 

Proof. Let us assume that there is an additional integral of motion J : P ^ 
R. By locahzing J around 7 if necessary, the Hamiltonian flow generated 
by J exists for all times s e R. Since these integrals are in involution, their 
flows commute: 

= *,$t (i,seR). 
Let T be the period of 7 and x e 7. We have 

$T*.(x) = ^^s^Ax) = ^s{x) (s e R). 

This implies that "ifg maps (within a given energy surface T^e) periodic orbits 
of the flow $t into periodic orbits of this flow. But 7 C is an isolated 
periodic orbit and therefore ^s{x) £ 7 for all s e R. This implies that H 
and J are functionally dependent on 7. □ 

Rerricirk 1 Note that single isolated periodic orbits need not form an ob- 
struction to the existence of an additional analytic integral of motion J : 
P — > R, independent of H in the sense defined in Sect. 4. If the periodic 
orbit 7 is hyperbolic, then J must be constant on its stable and unstable 
manifold, but still on a neighbourhood of 7 C Eg the singular set of J|s^ 
may only consist of 7 which is of measure zero. 

Easy examples for this are given by motion in a smooth potential V : 
R^ — > R which are rotationally symmetric. Then angular momentum is an 
independent analytic integral of the motion, but circular hyperbolic orbits 
may exist for some energies. 

In the present context the energy surface (£" > ^'th) for the two-centre 
problem contains exactly one bounded orbit 7. 7 is periodic and collides with 
the two centres. As commented in the Introduction, the two-centre problem 
is analytically integrable. 
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Let us return to the n-centre problem. By Thm. 12.8 from |Knj the set 
of bounded orbits Be on the energy level E > Eti^ consists in hyperbolic 
trajectories, and for n > 3 is locally homeomorphic to a product of a Cantor 
set and the interval. Moreover its Liouville measure vanishes. We now show 
that bE can form an obstruction to the existence of independent analytic 
integrals: 

Theorem 2 If the ajfine span of the (non-collinear) centres Si, . . . , s„ equals 
M.^, then for E > E' > E'th the n-centre problem does not admit a pair 
of independent analytic integrals /i,/2 : ^e 1R of motion on the energy 
surface T^e '■= H^^{E). 

Proof. Assume that Ji, J2 : ^ M are analytic integrals of motion. Then 
on the five-dimensional manifold S^; we consider the singular set 



with / := (/i,/2). 

By Prop, mthe periodic orbits within bs belong to the (closed) set S (as 
in Prop. E] we consider functions J : P — M, we use appropriate smooth 
extensions J of J/; : E^; M) . 

In jKn] it was shown that for energies E > E'th the periodic orbits are 
dense within bE- Thus bE C S, too. 

According to |vdDj (see also pjij where a sketch for more general suban- 
alytic sets appears) S admits an analytic simplicial decomposition which is 
locally finite and whose simplices are semianalytic. As bE <Z S is compact, 
there exist (disjoint) simplices Ai, . . . , with 



Now we assume /i,/2 to be indendent integrals, contrary to the statement 
of the theorem. Then K := maXjdim(Aj) < dim(E£;) = 5. As the set S is 
<I>j-invariant, we consider the transversal intersections 



of the simplices with the Poincare surface He C T,e (defined in Sect. 10 
of ^). We have K := maxidim(Ai) = K - 1 < dim(HE) = 4. The 



5 := {a; G I rank(D/(x)) < 2} 



m 




(22) 



A, := A, n He 



(23) 
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intersection Ae := Be (^'He of the bounded orbits with the Poincare surface 
has the form 

As = n A^, 

Ag C He being the (un-) stable manifolds, consisting of two-dimensional 
leaves which intersect transver sally. 

It is known from Hasselblatt [B.a\ that the Holder regularity of (the distri- 
butions of) A^ can be controlled by the so-called bunching constant. For the 
case of the n-centre problem Prop. 11.2 of |Knj controls the expansion resp. 
contraction rates of the Poincare map on A^, which differ from a constant 
times E by 0{E^). Thus applying [Ha] we have C^"*^ regularity of A^ for all 
large energies E > Eth- 

Now as i^T < 4, for each Aj in at least one of the intersections 

Af :=A,nA| 

must be of dimension < 1. By reversibility of the flow we assume w.l.o.g. 
that dim(A]~) < 1 and derive a contradiction. 

As A^; is a Cantor set (see Thm. 12.8 of Kn ), we assume w.l.o.g. that 
A]" contains a sequence {x^^^)i^fq of points x*^*^ G A^; converging to x G A]". 
By going to a subsequence, if necessary, we assume that the unit vectors 
t>(*) G TxHe with exp^(t*^*^f *^*^) = a;^*-* have a limit v := limj^oo 'J^*-*'' (which is 
in fact independent of the choice of Riemannian metric used to define the 
exponential map). 

In |Knj starting from the alphabet S := {1, . . . ,n} the space of symbol 
sequences 

X := {^G 5^ I Vi G Z : h+i ^ h] 
is equipped with the metric 

^^ai):= 5^2-1*1 -(l-^^/J, (fc,/GX). 

Then denoting the shift by 

a : X ^ X , cr(^)j := A;j+i (i G Z), 
there exists a Holder homeomorphism 

: X ^ As (24) 
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conjugating a and the restriction of the Poincare map to A^;, see Lemma 
12.2 of jEn|. 

The Poincare surface is the disjoint union 

n 

riE — [J , 

k.l = l 

the being open regions in the intersection of a five-dimensional affine 
space and J2e- Starting with 

VE{ko,h) ■.= WE{ko,h) ■.= n'S''\ 

in )Knj for • • • , ^o) admissible (that is ki ^ the nested subsets 

WE{k.m. . . . , A;o) := WE{k_^, ^o) n VE{WE{k_^, . . . , (25) 

resp. for (/cq, • • • , km) admissible 

l/s(A;o, . . . , := VE{k^, k^) n VfiyE^k^, • • • , A;™)), (26) 

were defined, using the Poincare map Ve- 

By going to subsequences we can assume that the points x*^*-* G corre- 
spond to symbol sequences e X which are related to the symbol sequence 
k := T^{x) G X of X by 

^ = k, (zGN;j>x(0) 

but /c^i) 7^ ^x(i)' where x(«) ^ -oo. 

By construction the vector v is tangent to the one-dimensional manifold 
I\\ at X. 

We now show the existence of a second sequence (?/'-*'')iGN of points y*^*^ G 
A^; converging to x, but with the following property: There exists an a G 
(0, 1) such that writing the points in the form 

l/« = exp,(s«w;(*)) 

with units vectors w*^*^ G TxH-e-, 

rf(w;('\span(t;)) > (27) 
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Namely as the number n of centres is > 4 (which follows from our as- 
sumption on the positions Si, . . . , s„) we find symbol sequences G X with 

but /JJ.) 7^ k^{i), and 

affine span ( Sfe Sfc^ , s^w ) = Ml (28) 

Then the y*-*-* := J-'e{L^^^) converge to 

lim y« = J^E f lim ij'A = J^E{k) = x. 

i—*oo \i— >oo / 

Next we consider the geometric situation at the (early) time More 
precisely we set 

^:=a^«(^) , ^«.= ^xW(^«) and f ^= (T^« (/(*)). 

Thus 

but h ^ ^ if ^ ko- 

By the conjugacy property of pH) the identities x := J-'E{k) = 
:= :Fe {t^) = Vf\x^^y) and y := Te (f ^) = Vf^ (y^) are true. 
All of these points are contained in the local stable manifold 

VE{ki, h, . . .) c VEik, h) c 

but at the same time in the following disjoint sets 

xeWE{hMM) , x^'^ eWE{ko\~ku~k2) , y'-'^ eWE{ll;\~ki,~k2). 
In 7^^^'^^ C Ti-E the minimal angle between vectors from (points in) 

WE{ko,h,k2)nVEih,~k2,~k3) to WECko\h,k2)r]VECki,k2,h) 
and vectors from 

WECko,~ki,h)r]VE{h,~k2M to WE{lf,ki,h)r]VEih,k2,h) 
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is bounded away from zero by some A > (see (I26|) and (|25|)). This follows 
from (j28|) . As there are only finitely many (at most n, to be more precise) 
choices for each of the indices ko, ki, k2, k^^ and Iq\ this bound is uniform. 
So the angle between the unit vectors w*-*^ and w^^^ defined by 

is bounded below by A. But by Prop. 11.2 of jKnj the contraction rates of 
these vectors w.r.t. the Poincare Ve differ at most by the order 0{1/E). 

By estimate (11.2) of |Knj we conclude that after iterating — x(0 > 
times these vectors still have an angle bounded away from zero by 

whereas their length is reduced at least by a factor (ci?)'^'-*^'. 

Thus by enlarging the energy E, we can choose an arbitrarily small a > 
in (j77|). 

But ()27p is incompatible with y^'^^ belonging to the one-dimensional submani- 
fold A^. 

As we can apply the same argument to all x G A^;, we have derived a 
contradiction to the finite covering assumption (|22p . □ 

The second author (I.A.T.) was supported by RFBR (grant 03-01-00403) 
and Max-Planck-Institute on Mathematics in Bonn. 
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